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The Glasgow reweighting method is evaluated for SU(2) lattice gauge theory at nonzero ^ and finite T. We 
establish that the 'overlap problem' of SU(3) measurements, in which the transition points determined from 
thermodynamic observables have an unphysical dependence on the value of fj, used to generate ensembles for 
reweighting, persists for SU(2). By combining the information from different lattice ensembles we alleviate 
sampling bias in the fugacity expansion, and identify the Lee Yang zeros associated with the transition to a high 
density phase that can plausibly be associated with diquark condensation. We also confirm the existence of a line 
of first order transitions above a critical point in the T — /j, phase plane previously predicted by effective chiral 
lagrangian calculations. 



1. Introduction 

Recent speculation over BCS type-instabilities 
in the Fermi surface at high density [|] |^ has led 
to resurgence of interest in two colour QCD. Lat- 
tice models with pseudoreal representations are 
attractive candidates for simulation at finite den- 
sity since for such theories the Dirac matrix is 
positive definite, which permits the use of exist- 
ing Monte Carlo techniques at n ^ 0. For SU(3) 
at finite density the lattice action is otherwise 
complex and so the probabilistic importance 
sampling step of such methods is therefore unde- 
fined. Reweighting methods have proved a useful 
means of addressing this issue, where the fi de- 
pendence of the grand canonical partition func- 
tion Z{fi) can be made semi-analytic in a fugac- 
ity expansion, as with the Glasgow method [Q. 
The complex action issue is thus avoided, by gen- 
erating a lattice ensemble in an accessible regime 
of the parameter space (eg. /i = for SU(3)). 
Naively, one might anticipate that the specific 
lattice ensemble used in the reweighting has little 
impact on numerical evaluations of the expansion. 
In fact quite the reverse is true. For the SU(3) 
model with dynamical quarks, even at interme- 
diate coupling where one might expect thermal 
fiuctuations to enhance the frequency of sampling 
physically relevant states, an ensemble generated 
at /X = reproduces a similar phase structure to 
the quenched model This has severe con- 



sequences since quenched finite density QCD is 
understood to be the zero fiavor limit (n — > 0) of 
a theory with equal numbers of quarks and con- 
jugate quarks. The lowest lying baryonic state in 
the model is thus the unphysical "baryonic pion" , 
formed from quark-conjugate quark pairs rather 
than the hghtest three quark state How- 
ever, this is of no consequence in two colour QCD 
1^ as the baryonic pion and baryon propagators 
are equivalent at fi ^ 0. 

Pseudoreal models are not entirely free from 
the effects of the reweighting overlap patholo- 
gies, however. Models with quarks in the adjoint 
representation (with n odd) suffer from the re- 
lated reweighting pathology: the sign problem. 
Since detM(^) is always real for two colour QCD, 
although reweighting is no longer mandatory it 
provides a useful opportunity to investigate the 
overlap issue (where the correct physics may be 
more easily extracted from lattice measurements 
by conventional means), and to quantify the sig- 
natures of the sampling numerical discrepancies 
therein. 

2. Symmetries of the SU(2) Lattice Action 

The n— fiavor symmetry of SU(2) QCD given 
by quarks in their fundamental representation 2, 
is SU{2n) rather than U{n)L x U{n)R as might 
be anticipated, as is demonstrated by a change of 
basis in the free field langrangian. 
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Table 1 

Global n— flavor symmetries of the two colour 
QCD action, in the contimmm and with staggered 
quarks for m, ^ 0. 
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where -0 is a 2 doublet. The inclusion of ex- 
plicit symmetry breaking terms in m and /i can be 
shown similarly to lead to the symmetry breaking 
patterns tabulated above p[pO|. 
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Naturally, the lattice model (with Kogut- 
Susskind fermions), follows a somewhat similar 
symmetry breaking scheme, having a manifest 
global U{2n) symmetry. In the continuum limit, 
for the choice of n = 1, this lattice action cor- 
responds to 8 physical flavors through the well- 
known doubling of fermionic mode s fllfl . 

At m, /X 7^ it is argued in |9| that since 
the number of the Goldstone modes differs for 
nonzero expectations of t/jV ^nd that a phase 
transition occurs at /i > ^ttItt corresponding to 
the point at which the number density of quarks 
becomes nonzero and U{n)y is spontaneously 
broken. It is then further argued with the Lan- 
dau free energy in [p^ , that since the number of 
Goldstone modes at high density is odd, that the 
transition to the free quark phase at finite tem- 
perature is necessarily first order. 

3. Glasgow Method 

For the Glasgow reweighting method, the 
dependence of the lattice action is made ana- 



lytic through the formulation of the fugacity ex- 
pansion, where z = exp(/i/T). This constitutes 
the characteristic fugacity polynomial which is 
formed from the propagator matrix P, defined 
through the fermion matrix M Where P 

is written in terms of the matrices which contain 
links between lattice sites in the spatial directions 
G, and forward and backward in the time direc- 
tion V and V'^ , 



2iM = 2im + G + Vef + We"'^ 
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detM = det(G + 2mi + y^e-^ + ye^) (7) 
= e"<="'"*^ det(P - e-") 
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with nlrit the lattice volume and ric the number 
of colours in the expansion. Since V is an overall 
factor of P the order of the expansion is reduced 
by exploited the unitary symmetry Zn^ defined 
by multiplying the timelinks in V by e^irij/nt^ 
where j is an integer. Since rit = 1/T, the Grand 
Canonical Partition function Z{ii) is thus given 
defined as an expansion in terms of the fugacity 
variable and the canonical partition functions Z„. 
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By reweighting this expansion an arbitrary nor- 
malisation to Z{^) is introduced, though this 
leaves the analytic determination of thermody- 
namic variables unaffected. 
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However, the reliability of the ensemble-averaging 
is strongly affected when the ratio of the ratio in 
eqn.(^|) differs greatly from one. The canonical 
partition functions can, in general, only be reli- 
ably determined for a small series of terms in the 
expansion, centered on the term of order n(/io). 
This effect and the reliability of the averaging 
elsewhere can be established for two-colour QCD 
by measuring the ratio of the ensemble-averaged 
expansion coefficients between two or more en- 
sembles generated at different values of fio- Hav- 
ing then identified n{^o) for several ensembles, 
our composite reweighting method then consists 
of rescaling the expansion coefficients from dif- 
ferent ensembles through these ratios (where the 
ensemble-averaging is effective). The bias in- 
troduced through reweighting the expansion can 
thus be systematically alleviated, and thermody- 
namic observables more reliably determined Jlif . 

3.1. Thermodynamic Observables 

The eigenvalues A„ of P naturally share the 
symmetries of V, most notably A — s- 1/A* re- 
lating P to P~^ up to a unitary transformation. 
Since SU(2) with quarks in the fundamental rep- 
resentation is pseudoreal it can also be shown that 
A — > A*. By rewriting the expansion in the vari- 
able y — z + l/z the order can be further reduced 
by a factor of two to reduce rounding errors in the 
numerical implementation p^ . The quark num- 
ber density n and its associated susceptibility Xn 
for this expansion is then both easily evaluated 
from the expansion, and in addition is also read- 
ily amenable to composite reweighting approach 
described above. 
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Similarly the zeros a„ of Z{fi) are readily identi- 
fied from the expansion, which as Lee and Yang 



showed with an Ising ferromagnetic system, cor- 
respond to a phase transition in the thermody- 
namic limit wherever a zero approaches the real 
axis in the complex-z plane |16|. 
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4. Results 



4.1. Intermediate Coupling 

We generated a total of seven ensembles at con- 
secutive values of fio ranging from /Iq = 0.3 — 1.1, 
at Pc = 1.5 both for 4^ and 6^4 lattice vol- 
umes. From these lattice ensembles we evaluated 
the Lee- Yang zeros, quark number density sus- 
ceptibility (xra), and {il^ip) using a conventional 
stochastic approach. 

For both our measurements at (3 = 1.5 and 
(i — 2.3, {4>i^) decreases gradually to zero over 
the range of values of /i we generated for m = 
0.05. However, since ('i/'V' (/^ = 0)) is consider- 
ably smaller in the latter measurement, plausibly 
U{1)a is spontaneously broken in the chiral limit 
in the former case for the volumes we used. It 
then follows that there should be a correspond- 
ing transition in the m ~ fi plane at jic ~ \Tn-K 
(as we argued in Sec 2), which we were able to 
identify from our Lee- Yang zeros measurements 
using composite reweighting. An unphysical /io 
dependence dominates our measurements prior to 
composite reweighting at (3 — 1.5 and is tabu- 
lated in Table 2, along with the convergence of 
our measurements after composite reweighting as 
we increase the number of included ensembles. 

For an ensemble generated at jjLo = He we be- 
lieve the sampling should be effective enough to 
circumvent the need for composite reweighting. 
There is evidence to support this with the en- 
semble we generated at = 0.3 ~ pc in Fig-1, 
which shows more evidence of a transition (where 
the zeros consistently approach the real axis) at 
/ic ~ ^tItt than the other ensembles we gener- 
ated at (3 — 1.5. However, since we are unable 
to accurately quantify which values of HcifJ-o) are 
the more valid from our jacknife error estimates 
of the Lee- Yang zeros, and the unphysical /io de- 
pendence of our measurements persists for ensem- 



bles generated at values of arbitrarily close to 
fic, we found it is more effective to sample the 
expansion coefficients accurately by generating a 
covering series of ensembles with our composite 
reweighting method. 

In varying the lattice volume V and /3, we can 
confirm that this unphysical fio dependence in our 
measurements behaves as we would expect of a 
reweighting overlap (sign) problem. The expec- 
tation of the sign of the Monte Carlo measure 
(which is treated as an observable for reweighting 
in the Potts model), shows a /? and V dependence 
of the form, 



{sgn) = — - exp(-/3yA/) 
^1 
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where Z| | is the partition function of the ensemble 
modified to exclude the sign problem amenable 
to a Monte Carlo approach, and A/ the differ- 
ence in free energy densities between ensembles 
0. This eflFect is seen Tables 2 and 3 where 
the imaginary part of the zeros nearest the real 
axis evaluated from the ensemble at generated at 
Mo Mc is comparatively smaller (and therefore 
more convincing) as V is increased. Similarly, this 
Lee- Yang zeros effect becomes more pronounced 
as we increase f3. 

We are able to determine the range of values 
of /io over which to generate ensembles for our 
effective sampling strategy, from the jacknife er- 
ror estimates for the ensemble-averaged expan- 
sion coefficients, which give us n(/io)- The largest 
of the coefficients C2„^n3 (related to the canonical 
partition function for the filled lattice) is of or- 
der one for fio = 1.2, and the lattice therefore 
saturated. Our quark number density suscepti- 
bility measurements (Xn(M)) become singular at 
fic as we include more ensembles in the compos- 
ite reweighting across this range, indicating that 
the transition at /ic ~ for f3 = 1.5 is first 

order. Fig 4. We are also able to identified a sec- 
ond smaller peak in these measurements which 
corresponds to the point at which the expecta- 
tion of the diquark falls off in existing conden- 
sate measurements. As saturation is approached 
at /X = 1.2 the diquark condensate thus evapo- 
rates in a less well determined transition driven 
by Fermi statistics |lll[|9l. 
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Figure 1. Lee Yang zeros evaluated in the com- 
plex fi plane {ijn = Tin an) for a 6^ A lattice at 
(3 = 1.5 from an ensemble generated at jig ~ 0.3. 
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Figure 2. Lee Yang zeros evaluated in the com- 
plex /i plane for a 6'^4 lattice at /? = 1.5 from an 
ensemble generated at /io=0.5. 
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Figure 3. Lee Yang zeros evaluated in the com- 
plex /i plane for a 6'^4 lattice at (3 — 1.5 from an 
ensemble generated at fio—0.7. 
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Table 2 

Lee Yang zero with the smallest imaginary part 
evaluated in the complcx/i plane {rjn = Tin an) 
for two lattice volumes at /3 = 1.5. Dependence 
on value of /Xo used to generate ensembles for the 
Glasgow reweighting method (upper) , and depen- 
dence on the number of ensembles includc;d in the 
new composite reweighting scheme (lower). 




H„ = {0.5, 0.7', 0.8, 0.§, 1.0, 1.^, 1.2} 
H„ = { 0.3, 0.5, 0.7, 0.8, 0.9, 1 .0, 1 .1 , 1 .2, 1 .3 ) 
3.5 n„ = { 0.1 , 0.3, 0.5, 0.7, 0.8, 0.9, 1 .0, 1 .1 , 1 .2, 1 .3, 1 .4 J 

t 



2.5 
2 

1.5 
1 

0.5 




/hi 



0.2 



0.4 



0.6 



0.8 



1.6 



Figure 4. Quark number density susceptibility 
(Xn(A*)) for a 4^ lattice at /3 = 1.5, with an 
increasing number of ensembles included in the 
composite reweighting (upper 1, 3, 5) (lower 7, 
9, 11). A prominent peak develops as the num- 
ber of composite reweighted splines is increased, 
indicating a first order transition. 



4.2. Weak Coupling 

There is a niarlced difference between the un- 
physical fio dependence of the Lee- Yang zeros 
associated with fic evaluated before composite 
reweighting at (3 = 1.5 and (3 — 2.3 between Ta- 
bles 2 and 3. In Table 2 there is some indication of 
competition between the two separate transition 
points during rootfinding we have identified above 
as the value of fio is varied, which is now entirely 
absent in Table 3. From this we can conclude that 
there is no indication of a second transition point 
at /3 = 2.3, and that a transition can be read- 
ily identified at Hc ~ 0.8. Where the reweight- 
ing ensemble is generated at /io — fic there is 
also good agreement between these zeros mea- 
surements and those from composite reweighting. 
Although again we have not quantified how small 
the difference between fio and fic must be for the 
reweighting to be effective, and have relied in- 
stead on composite reweighting. 

Our composite reweighting Lee- Yang zeros and 
quark number density susceptibility measure- 
ments indicate again (where the zero closest the 
real z-axis goes to zero as the volume is increased 
and where (XnifJ-)) becomes singular) that the 
transition is first order. The only context in which 
a first order transition is predicted in the effec- 
tive chiral lagrangian approach is with a transi- 
tion from the diquark to the symmetric phase for 
/ic > ^Wtt, and we therefore confirm the existence 
of such a transition line in the T — ji plane. 

5. Conclusions 

Despite expectations, generating an ensemble 
for the reweighting method with a value of /io ar- 
bitrarily close to fic still leads to an overlap prob- 
lem in a model with a pseudoreal representation 
evaluated at (3c- In fact, for an exploratory study 
(in which fic is unknown), the fugacity expansion 
coefficients are more effectively sampled through 
the combination of terms from a covering series 
of ensembles. Even with the real Monte Carlo 
measure of two colour QCD the overlap prob- 
lem is still pathological for the Glasgow reweight- 
ing method. With new multi-parameter (/?, fj.) 
reweighting approaches to SU(3) |Q we would 
therefore expect there to be similar sampling bias 
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Figure 5. Lee Yang zeros evaluated in the com- 
plex fi plane for a 6'^4 lattice at (3 = 2.3 from an 
ensemble generated at fio~0-7. 
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Figure 6. Lee Yang zeros evaluated in the com- 
plex fi plane for a 6'^4 lattice at (3 = 2.3 from an 
ensemble generated at fio— 0-8 ■ 
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Figure 7. Lee Yang zeros evaluated in the com- 
plex fi plane for a 6'^4 lattice at (3 — 2.3 from an 
ensemble generated at /io=0.9. 
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Table 3 

Lee Yang zero with the smallest imaginary part 
evaluated in the complex/i plane (ry„ = Tin an) 
for two lattice volumes at /3 = 2.3. Dependence 
on value of fio used to generate ensembles for the 
Glasgow reweighting method (upper), and depen- 
dence on the number of ensembles included in the 
new composite reweighting scheme (lower). 



in the coefficients across the full fugacity expan- 
sion. In order to successfully investigate the pos- 
sibility of there being similar tricritical behavior 
in the /i — T phase plane of SU(3) with the 
Glasgow method, this issue of effective sampling 
across the full range of the polynomial this over- 
lap problem should then be addressed. As we 
have, seen substantially different transition points 
are determined from reweighting measurements 
at intermediate coupling from ensembles gener- 
ated with Ho arbitrarily close to fic- 

Lee Yang zero analysis allows the simple iden- 
tification of a first order transitions in jj, with 
lattice measurements on comparatively small vol- 
umes. To extend the rigor of this approach and 
determine the critical exponents of the transition 
at fic, however, it will be necessary to increase 
the lattice size as the zeros scaling at /3 = 2.3 is 
not without finite volume effects. It will also be 
interesting to repeat this volume scaling analysis 
for (3 < Pc and to compare the measured crit- 
ical exponents for the diquark phase transition 
with those predicted with the chiral langragian 
approach |l^, believed to be second order. 

Thanks to M. Alford for useful discussions. 
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